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The Durham/UKST Galaxy Redshift Survey - IV. 
Redshift Space Distortions in the 2-Point Correlation 
Function. 
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ABSTRACT 

We have investigated the redshift space distortions in the opticaUy selected 
Durham/UKST Galaxy Redshift Survey using the 2-point galaxy correlation func- 
tion perpendicular and parallel to the observer's line of sight, ^(cr, tt). We present 
results for the real space 2-point correlation function, ^(r), by inverting the optimally 
estimated projected correlation function, which is obtained by integration of ^(cr, tt), 
and find good agreement with other real space estimates. On small, non-linear scales 
we observe an elongation of the constant ^(cr, tt) contours in the line of sight direction. 
This is due to the galaxy velocity dispersion and is the common "Finger of God" ef- 
fect seen in redshift surveys. Our result for the one-dimensional pairwise rms velocity 
dispersion is < w^ >i/2— 416 ± 36kms~^ which is consistent with those from recent 
redshift surveys and canonical values, but inconsistent with SCDM or LCDM models. 
On larger, linear scales we observe a compression of the f(cr, tt) contours in the line 
of sight direction. This is due to the infall of galaxies into overdense regions and the 
Durham/UKST data favours a value of {H.^'^ /b)^0.5, where O is the mean mass den- 
sity of the Universe and b is the linear bias factor which relates the galaxy and mass 
distributions. Comparison with other optical estimates yield consistent results, with 
the conclusion that the data does not favour an unbiased critical-density universe. 

Key words: galaxies: clusters - galaxies: general - cosmology: observations - large- 
scale structure of Universe. 



1 INTRODUCTION 

The use of redshifts as distance estimates is commonplace 
in surveys which map the galaxy distribution of the Uni- 
verse. While redshifts are both quick and easy to acquire, 
they do not generally reflect the true distance to the galaxy 
because of the galaxy's own peculiar velocity with respect 
to the Hubble flow. Therefore, the observed clustering pat- 
tern will be imprinted with the galaxy peculiar velocity field 
along the observer's line of sight direction. While at first 
this would appear to be highly problematic for the use and 
understanding of the measured clustering statistics, these 
so-called redshift space distortions can be used to estimate 
some important cosmological parameters describing the dy- 
namics of the Universe (e.g. Peebles 1980; Kaiser 1987). 

We will be investigating these redshift space distortions 
using the spatial 2-point correlation function, ^, and denote 
a real space separation by r and a redshift space separa- 
tion by s. Non-linear and linear efi'ects will be competing 
on all scales but for the most part we take the simplifying 



assumption that non-linear effects dominate on small scales 
(< 10/i~^Mpc), whereas linear effects dominate on larger 
scales (> 10/i~^Mpc). Although this is a slightly naive ap- 
proach we use N-body simulations and mock catalogues to 
guide us in determining the accuracy of the modelling in- 
volved. 

The initial clustering results, redshift maps, etc. of the 
Durham/UKST Galaxy Reds hift Survey were sum marized 
in the first paper of this series (Flatcliffe et al. 1996a). In this 
paper we present a detailed analysis of the redshift space dis- 
tortions in the 2-point correlation function as estimated from 
this optically selected survey. We briefly describe our survey 
in Section Section gives a qualitative description of the 
effects of the redshift space distortions on the 2-point corre- 
lation function. In Section we describe and test methods 
of obtaining the real space correlation function from red- 
shift space estimates. In Section |5| we describe non-linear 
effects and estimate the one-dimensional rms galaxy pair- 
wise velocity disperison. Linear infall effects are described 



in Section M and an estimate of Q /b is obtained. The re- 
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suits of this analysis are discussed and compared with other 
redshift surveys and models of structure formation in Sec- 
tion M. Finally, we summarize our conclusions in Section H. 



2 THE DURHAM/UKST GALAXY REDSHIFT 
SURVEY 

The Durham/UKST Galaxy Redshift Survey was con- 
structed usin g the FLAIR fibre optic system ( Parker & 
Watson I995I ) on the 1.2m UK Schmidt Telescope at Sid- 
ing Spring, Australia. This survey uses the astrometry 
and photometry from the Edinburgh/Durham Southern 
Galaxy Catalogue (EDSGC; Collins, Heydon-Dumbleton 
& MacGiUivray 1988; Colhns, Nichol & Lumsden 1992) 
and was completed in 1995 after a 3-yr observing pro- 
gramme. The survey itself covers a ~20° x 75° area cen- 
tered on the South Galactic Pole (60 UKST plates) and 
is sparse sampled at a rate of one in three of the galax- 
ies to bj ~ 17 mag. The resulting survey contains ^2500 
redshifts, probes to a depth greater than SOO/i'^Mpc, with 
a median depth of ~150/i^^Mpc, and surveys a volume of 
space ~4 x Vf[h~^Mpcf. 

The survey is >75 per cent complete to the nominal 
magnitude limit of bj = 17.0 mag. This incompleteness 
was mainly caused by poor observing conditions, intrinsi- 
cally low throughput fibres and other various observational 
effects. In a comparison with ~150 pubhshed galaxy veloci- 
ties (Peterson et al. 1986; Fairall & Jones 1988; Metcalfe et 
al. 1989; da Costa et al. 1991) our measured redshifts had 
negligible offset and were accurate to ±150 kms~^. The scat- 
ter in the EDSGC magnitudes has be en estimated at ±0.22 
mags (Metcalfe, Fong & Shanks 1995) for a sample of ~100 



galaxies. This scatter has been confirmed by a preliminary 
analysis of a larger sample of high quality CCD photometry. 
All of these observational details are discussed further in a 
forthcoming data paper (Ratcliffe et al., in preparation). 



3 REDSHIFT SPACE DISTORTIONS 

We investigate the effects of redshift space distortions by 
estimating the spatial 2-point correlation function, ^, as a 
function of the two variables a and vr. These are the separa- 
tions perpendicular {a) and parallel (vr) to the line of sight. 
The specific definitions of a and tt we use were given in Rat- 
cliffe et al. (1996c), hereafter Paper III. However, we found 
that our results do not depend significantly on their exact 
form. 

Galaxy peculiar velocities cause an otherwise isotropic 
real space correlation function to become anisotropic when 
observed in redshift space. The degree of anisotropy mea- 
sures the low-order moments of the peculiar velocity distri- 
bution function (e.g. Peebles 1980). On small, non-linear 
scales (< 10/i~^Mpc) the velocity dispersion of galaxies 
in virialised regions dominates the anisotropy, causing an 
elongation in the contours of constant ^ along the line of 
sight direction. This is the well-known "Finger of God" ef- 
fect seen in redshift surveys and allows one to estimate the 
one-dimensional pairwise rms velocity dispersion of galaxies, 
< 'Up >^'^ (e.g. Davis & Peebles 1983). On larger, more lin- 
ear scales (> 10/i~^Mpc) the infall of galaxies into overdense 



regions dominates, causing a compression of the ^ contours 
in the line of sight direction (Kaiser 1987). This effect al- 
lows a measurement of f2°'^/&. While our division of the 
non-linear and linear regimes is slightly naive, we will be 
guided by CDM mock catalogues to determine the accuracy 
of the modelling in the different regions. 



3.1 The Redshift Space 2-Point Correlation 
Function, ^(f, tt) 

We have estimated ^((t, tt) using the optimal methods de- 
termined empirically in Paper III. Briefly, this involves dis- 
tributing a random and homogeneous catalogue with the 
same angular and radial selection functions as the original 
survey. One then cross correlates data-data, data-random 
and random-random pairs, binning them as a function of 
separation (in this case a and tt). As a result of testing 
Monte Carlo mock catalogues drawn from CDM A''-body 
simulations we found that the estimator and weighting com- 
bination which most accurately traced the actual ^ and also 
produced the minimum variance in ^ was the estimator of 
Hamihon (1993) and the weighting of Efstathiou (1988). All 
of these techniques and the biases in them are discussed at 
length in Paper III and we refer the interested reader there. 
In Figs, hi and EI we show contour plots of constant ^ as 
a function of a and tt. On both of these figures we adopt the 
following conventions: solid lines denote ^ > 1 with A^ = 
1; short-dashed lines denote < ^ < 1 with A^ = 0.1; 
and long-dashed lines denote ^ < with A^ = 0.1. For 
reference, the contours ^ = 1 and are in thick bold and the 
regularly spaced thin bold lines show an isotropic correlation 
function for comparison. When the plotted scale is log- log 
the binning size is 0.2dex, when it is linear-linear the binning 
is 1.0/i~'^Mpc. In both cases no formal smoothing has been 
applied. 



3.2 Results from the CDM models 

We have calculated the ^(cr, vr) results from two different 
CDM models (Efstathiou et al. 1985; Gaztafiaga & Baugh 
1995; Eke et al. 1996): standard CDM with Qh = 0.5, b = 1.6 
(SCDM); and CDM with Qh = 0.2, 6 = 1 and a cosmolog- 
ical constant (A — 0.8) to ensure a spatially flat cosmology 
(LCDM). In Fig. Wa) we show the average ^(ct, tt) calculated 
directly from the 9 SCDM A^-body simulations which were 
available to us, while Fig. [y(b) shows the average from the 5 
available LCDM simulations. Given that these simulations 
are fully volume limited with a well-defined mean density 
there are no problems with estimator/weighting schemes 
and we estimate ^ using the methods described in Paper 
III for this case. For simplicity, we also use the distant ob- 
server approximation (i.e. we imagine placing the A^-body 
cube at a large distance from the observer) and therefore 
the line of sight direction can be assumed to be one specific 
directi on. We choose this to be the z-direction and hence 
a = y'x'^ +y^, TV = z and the bins are cylindrical shells. 

Fig. Ka) shows that the small scale velocity disper- 
sion in the SCDM model dominates the whole plot, elon- 
gating the contours even on large scales (tt > Wh~ Mpc). 
The ^ = 1 contour cuts the a axis between 4.5-5.0/i~^Mpc, 
which agrees well with the real space correlation length of 
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a (h^'Mpc) 




a (h^'Mpc) 

Figure 1. The results of ^{a, w) calculated directly from Af-body 
simulations of two different CDM models. Fig. (a) shows the av- 
erage of the 9 available SCDM simulations, while Fig. (b) shows 
the average of the 5 available LCDM simulations. The diffe rent 
types of contours are described in detail in the text of Section p.l[ 
It is clear that the elongation of the contours in the line of sight 
(vr) direction, caused by the non-linear galaxy velocity dispersion, 
dominates the SCDM model. A similar effect is seen in the LCDM 
model, although at a lower level. 



plot, albeit less pronounced than the SCDM model. Indeed, 
there is possible evidence of a flattening in the n direction 
on ~20/i~^Mpc scales. The ^ = 1 contour cuts the a axis be- 
tween 6.0-6.5/i~^Mpc, which again agrees well with the real 
space correlation length of ~6.0/i~^Mpc estimated in Paper 
III. Again, there is a strong signal for the non-linear results 
to model but this time it appears possible that a sensible 
linear result could be obtained from the LCDM model (see 
Section 6). 

3.3 Results from the Durham/UKST Survey 

In Fig. y we show the results of ^(ct, vr) from the 
Durham/UKST survey as calculated using the optimal 
methods and magnitude limits described in Ratcliffe et al. 
(1996b) and Paper III. As mentioned in Section Bjl this in- 
cludes the estimator of Hamilton (1993) and the weighting 
of Efstathiou (1988). It is obvious that the noise levels for 
the data are significantly higher than those of the iV-body 
simulations. Given the differences in signal between the N- 
body simulations and the redshift survey this was to be ex- 
pected. As an aside, we note that the CDM mock catalogues 
used later in this paper have similar noise levels to the data 
in plots like these. Fig. ti shows that the non-linear veloc- 
ity dispersion does not dominate the whole plot (unlike the 
CDM models of Fig. hi) and there is no significant elongation 
beyond 7r~5-10/i~^Mpc. Also, on larger scales than these a 
visible compression of the contours in the n direction in seen. 
Therefore, we are confident that believable non-linear and 
linear results can be obtained from the Durham/UKST sur- 
vey using the approach attempted here. 



4 THE REAL SPACE 2-POINT 
CORRELATION FUNCTION 

We define the projected 2-point correlation function, ui„ (a) , 
by (e.g. Peebles 1980) 



Wv{o-) 



^{a,Tv)dTT, 



= 2 



^{y/a^ + n^)dn, 



(1) 
(2) 



where ^{\/a^ + vr^) is the real space correlation function. We 
calculate Wv{a-) using 



Wv{cr) = 2 / C(a, ■7r)d7r, 



(3) 



where the noise in ^(ct, tv) at very large scales makes us trun- 
cate the integral at some finite limit ncut- In practice we use 
a TTcut of 30/i~'^Mpc and our results are insensitive to raising 
this value. 



~5.0/i~^Mpc estimated in Paper III. While there is a strong 
signal to be modelled for the non-linear results, it is very 
doubtful that any useful information will be obtained for 
the linear results of the SCDM model without a more so- 
phisticated approach than is attempted here. 

Fig. p|(b) shows that the small scale velocity dispersion 
in the LCDM model is also the dominant feature on this 
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mri^k 




a (h^'Mpc) 




Figure 2. The results of {(c, it) estimated from the 
Durham/UKST survey using the optimal methods involving the 
estimator of Hamilton (1993) and the weighting of Efstathiou 
(1988). Fig. (a) plots the results on a log-log scale to empha- 
size the large scale features, while Fig. (b) plots the results on 
a linear-linear scale to emphasize the small scale features. The 
differen t ty pes of contours are described in detail in the text of 
Section |3.l| . One immediately sees that the noise levels in these 
diagrams are significantly higher than those of the A^-body sim- 
ulations shown in Fig. hi. These contour plots also show that the 
non-linear velocity dispersion elongates the contours in the line of 
sight (tt) direction out to 5-10/i~^Mpc. However, on larger scales 
there is evidence for a measurable flattening of the contours in 
the TT direction. 
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Figure 3. Testing the methods of estimating ^(r) using the Abel 
inversion of lu^ (ct) from (a) the SCDM and (b) the LCDM mock 
catalogues. We show a comparison of the actual ^(r) calculated 
directly from A^-body simulations (solid line) with the average 
5(r) obtained from the mock catalogues (solid symbols). These 
points came from taking each catalogue's optimally estimated 
5((T, vr), performing the integral in equation H to give a set of 
»i)^((T)'s, then carrying out the summations in equation a to give 
a set of f (r)'s and finally averaging to obtain the results. The error 
bars show the Icr standard deviation on this mean assuming that 
each mock catalogue provides an independent estimate of ^(r). 
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As we have discussed in Section pi any catalogue which 
uses redshifts to estimate distances suffers from the effects of 
galaxy peculiar velocities. Therefore, only the redshift space 
2-point correlation function, ^(s), is directly observable from 
a redshift survey. However, in terms of measuring and com- 
paring the clustering properties of the galaxy distribution 
with models of structure formation, the object of fundamen- 
tal interest is the real space 2-point correlation function, 
^(r). While (,{r) is not directly observable from a redshift 
survey we will see that it is possible to estimate it in an 
indirect manner. In this context we follow on from results 
presented in Paper III. 

4.1 The Projected 2-Point Correlation Function 

Our methods of estimating Wv(a-) were tested in Paper III 
where the results from the Durham/UKST survey were first 
presented (using the optimal methods mentioned in Sec- 
tion 3.3 for calculating ^). We then went on to model Wv{o') 
using a power law for ^(r) — (ro/r)^ and estimated the 
real space correlation length, ro = 5.1 ± 0.3, and slope, 
7 = 1.6 ± 0.1. In this section, rather than assume a power 
law form, we consider a direct inversion of equation to give 

4.2 Abel Inversion 

Equation H can be mathematically inverted using the gener- 
alized Abel equation to give (e.g. Lilje & Efstathiou 1988) 



1000 



ar) 



d[wy{a)] da 



da ^^ 



(4) 



Saunders, Rowan- Robinson & Lawrence (1992) consider the 
case when the data is logarithmically binned and Wv (a) then 
takes the form of a series of step functions, ■Wv{ai) — Wi, 
with logarithmic spacing centered on ai . They then linearly 
interpolate between each Wv point to get around singularities 
in the integral. In this case the d[wv{a)]/da factor simplifies 
to a constant value for each pair of ai spacings. The rest of 
the integral can be evaluated to give ^(r) at r = cr,; 



CC'^O 



TV ^-^ 

j>i 



Wj + l 



o-j+i 



ln|^^li±^^Si).(5) 



a-j + y/^ 



We test this method using mock catalogues drawn from 
the SCDM and LCDM A^-body simulations. These mock 
catalogues have the same angular/radial selection functions 
and completeness rates as the Durham/UKST survey and 
their construction was described in detail in Paper III. We 
selected the mock catalogues to sample independent volumes 
of the simulations, giving a total of 18 and 15 to analyse from 
the SCDM and LCDM simulations, respectively. The results 
of Wv{a) from these two sets of mock catalogues were pre- 
sented in Paper III and Fig. H shows the results of applying 
equation H to this Wv (a) data. The solid line shows the ac- 
tual ^(r) calculated directly from the A'^-body simulations 
and therefore denotes the answer we are trying to obtain. 
The solid symbols show the average ^(r) obtained from the 
mock catalogues by taking each catalogue's optimally esti- 
mated £,{a, n), performing the integral in equation ^ to give 
a set of w„(cr)'s, then carrying out the summations in equa- 
tion U to give a set of ^(r)'s and finally averaging. The error 



— (4.8/r)i-6 

• Durham/UKST 




(h-iMpc) 



Figure 4. Estimates of ^(r) from the Durham/UKST survey ob- 
tained by Abel inversion of Wv(a). The points are calculated by 
taking the optimally estimated ^(a, tt) from Fig. bl evaluating the 
integral in equation ti to give Wv{(t) and then inverting this with 
equation H to produce ^(r). The solid line shows the best fitting 
power law to the estimated §(r) in the indicated range. 



bars show the la standard deviation on this mean assuming 
that each mock catalogue provides an independent estimate 
ofe(r). 

The SCDM resuhs in Fig. Q{a.) show that the mock 
catalogues reproduce the actual ^(r) quite well out to 
~25/i~^Mpc scales (with some degree of scatter about the 
answer). On scales larger than this ^(r) is zero and there- 
fore there is no signal to invert, other than the noise in 
C(cr,7r). The LCDM resuhs in Fig. ||(b) show that the 
mock catalogues reproduce the actual ^(r) very well out to 
~40/i~^Mpc (apart from the point at ~7/!.~^Mpc which is 
caused by a couple of negative ^(r) points produced in the 
inversion process). Overall, we see that this method of in- 
version does reproduce the actual ^(r) and its features quite 
well, albeit with considerable scatter in places. 

This method of estimating ^(r) is applied to the 
Durham/UKST survey £,{a, vr) data of Fig. H. The w^ (a) re- 
sults were shown in Paper III and are used in equation M to 
produce the ^(r) shown in Fig. W. The error bars on this plot 
are obtained by splitting the Durham/UKST survey into 4 
quadrants, applying the method to each quadrant and then 
estimating a standard deviation of these results. These error 
bars are of a similar size to those estimated on an individual 
LCDM mock catalogue. The solid line shows the minimum 
X^ fit of a power law ^(r) = (ro/r)^ in the 0.1-20/i"^Mpc 
region. The best fit parameters were ro — 4.8 ± 0.5/i~'^Mpc 
and 7 = 1.6 ± 0.3, which gave a x^ of ~8 for 15 degrees of 
freedom. Note that the negative point at ~4.5/i~^Mpc was 
removed from this x^ statistic to avoid biasing the fit. Errors 
on these parameters come from the appropriate A^^ contour 
about the minimum. However, given the correlated nature 
of these points we anticipate that our quoted errors are more 
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than likely an underestimate. This should be adequate for 
the simple comparison done here. 



4.3 Richardson-Lucy Inversion 

Richarson-Lucy inversion has recently been made popular 
in large-scale structure applications by Baugh & Efstathiou 
(1993) and here we apply it to estimate ^(r) from Wv(o). 
Equation eI can be re- written as 



■Wv{o) = / i{r) 



2r 



\/r'^ — a^ 



dr. 



(6) 



by changing the variable of integration from n to r. If we 
then define 



K(a,r) = 



2r 



our integral becomes 



uii,(cr) = / ^{r)K{a,r)dr, 



for < r < cr, (7) 

for a < r < oo, (8) 

(9) 



which is in the form of a Fredholm integral equation of the 
first kind 



x) ^ tp{t)P{x\t)dt, 



(10) 



where (/"(k) is the known (or observed) function, ■(/'(i) is the 
unknown function and P{x\t) is the kernel of the integral 
equation. Therefore, equations can be numerically inverted 
using the method developed independently by Richardson 
(1972) and Lucy (1974). Basically, one guesses an initial 
form for V(i)i predicts (j){x) and uses this to update the 
original ip{t). This iteration continues and after n iterates 



(11) 



(12) 



0"(a;) = / i:"{t)P{x\t)dt, 

J a 

with the next iterate of ^^{t) given by 



r''\t) = r{t) 



I 






P{x\t)dx 



£p{x\t)di 



where (pi^) is the actual observed function. To apply this 
method to the logarithmically binned Wjj{a) data the inte- 
grals in equations |ll| and |l^ are approximated by the fol- 
lowing summations 



AT 



<i'^j) = > ^e{u)Kia„n)uAlnr, (13) 

e^\n) = Cin) ^';::^'^;] '' [' , (14) 

where M is the number of Wv (cr) bins and A'^ — A4/2 is the 
number of ^(r) bins. Our spacing in cr is Alger = 0.1 and 
hence our spacing in r is Algr = 0.2. Obviously one cannot 
get back more data points than are put in and A'' < M. 
In general, our choice of A'' = A//2 should assure a fairly 
smooth answer. 

There are two points worth noting about Richardson- 
Lucy algorithms. Firstly, there is no constraint on how many 



iterations are required for convergence to a stable answer. 
Therefore, there is no specific rule to know when to stop 
iterating. Too many iterations will cause convergence to the 
small scale noise features in the data, while too few implies 
that the results have not yet converged to the larger scale 
signal. Experience with Richardson-Lucy techniques shows 
that ~10 iterations are generally required (e.g. Lucy 1994). 
Secondly, this method assumes that the function tp{t) > 0. 
This is not always the case for our function ^(r). However, 
this is not too worrying as ^(r) is only likely to go negative on 
large scales when it is nearly zero, this is where our inversion 
process will be least believable anyway. Also, other authors 
(e.g. Baugh & Efstathiou 1993; Baugh 1996) have applied 
similar inversion techniques to the angular correlation func- 
tion to estimate the power spectrum (always positive) and 
real space correlation function (negative tail) and find very 
consistent results. 

We again test these methods using the SCDM and 
LCDM mock catalogues. Fig. H shows the results of apply- 
ing equations [l3| and p^ to the w„(cr) data. Again, the solid 
line shows the actual ^(r) calculated directly from the N- 
body simulations and the solid symbols show the ^(r) ob- 
tained from the mock catalogues by individually inverting 
Wv{cr) with the iteration process and averaging at the end. 
The circles show the results on stopping after 10 iterations, 
while the triangles show the results after 20 iterations. The 
results were found to be independent of the initial form of 
^(r), therefore a simple power law was used. The error bars 
show the lc7 standard deviation on this mean assuming that 
each mock catalogue provides an independent estimate of 
C(r). 

The SCDM results in Fig. |(a) show that the mock 
catalogues reproduce the actual ^(r) very well out to 
~30/i~^Mpc scales. On larger scales noise in the data dom- 
inates and this is the limit where we believe our inversion 
technique. It appears that the results do not change very 
much with further iteration. The LCDM results in Fig. Hb) 
show a similar form to those of the SCDM mock catalogues 
and reproduce the actual ^(r) very well out to ~40/i~^Mpc 
scales. Again, on larger scales noise in the data dominates 
and this is our believable inversion limit. Further iteration 
does not damage the stability of the inversion although the 
results do become less smooth. We note that this inversion 
technique does not appear to do well at the end-points of 
the integration range, this is an artifact of the finite inte- 
gration range (which assumes that the signal is zero outside 
of this range) and the large bins used in the summations. 
Therefore, our results could be improved by interpolating 
between the Wv{a) points and then inverting this function 
given that we could calculate it at many points. Given the 
uncertainties involved in ^(cr, tt) this is not attempted here. 
Overall, this method of inversion reproduces the actual ^(r) 
and its features very well and is arguably smoother than 
the Abel inversion. However, this is offset by only having an 
estimate at half the number of r points. 

This method of estimating f(r) is applied to the 
Durham/UKST survey £,{a, vr) data of Fig. 3. The Wv{a) re- 
sults were shown in Paper HI and are used in equations [13 
and [ij to produce the ^(r) shown in Fig. |6[ The error bars on 
this plot are obtained by splitting the Durham/UKST sur- 
vey into 4 quadrants, applying the method to each quadrant 
and then estimating a standard deviation of these results. 
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Figure 5. Testing the methods of estimating ^(r) using 
Richardson-Lucy inversion of ui„(fT) from (a) the SCDM and (b) 
the LCDM mock catalogues. We show a comparison of the actual 
^(r) calculated directly from A^-body simulations (solid line) with 
the average ^(r) obtained from the mock catalogues (solid sym- 
bols). These points came from taking each catalogue's optimally 
estimated ^{a,n), performing the integral in equation H to give 
a set of Wv(crYs. then carrying out the iteration between equa- 
tions hd and h4^ to give a set of ^{r)'s and finally averaging to 
obtain the results. The circles show the results on stopping after 
10 iterations, while the triangles show the same results after 20 
iterations. The error bars show the la standard deviation on this 
mean assuming that each mock catalogue provides an indepen- 
dent estimate of 5(r). 
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Figure 6. Estimates of ^(r) from the Durham/UKST survey ob- 
tained by Richardson-Lucy inversion of Wv{cr). The solid points 
are calculated by taking the optimally estimated ^(a, n) from 
Fig. y, evaluating the integral in equation M to give Wv{u) and 
then inverting using the iteration process involving equations [13 
and h4 to produce ^(r). The circles show the results on stopping 
after 10 iterations, while the triangles show the results after 20 
iterations. The solid line shows the best fitting power law to the 
estimated ^(r) in the indicated range. 



These error bars are of a similar size to those estimated on an 
individual LCDM mock catalogue. It can be seen that fur- 
ther iteration does not improve the results, in fact it makes 
them slightly less smooth. The solid line shows the minimum 
X^ fit of a power law ^(r) = (ro/r)^ in the 0,l-20/i"^Mpc 
region. The best fit parameters were ro = 4,6 ± 0.6/i~^Mpc 
and 7 = 1,6 ± 0,1, which gave a x^ of ~2 for 8 degrees of 
freedom. Note that the point at ~4,5/i~^Mpc (which was 
negative for the Abel inversion) is still low. However, it does 
not bias the fit because the larger binning in this method 
has essentially averaged over it. Errors on these parameters 
come from the appropriate Ax^ contour about the mini- 
mum. We make a similar comment to before about the cor- 
related nature of these points implying an underestimated 
error on these parameters. Again, this should be adequate 
for the simple comparison done here. 



4.4 Discussion 

Before discussing the results in detail we make a few com- 
ments about the inversion methods themselves. Firstly, 
given that the Abel inversion is a point by point process, 
the noise in w„(ct) is also inverted. This is similar to what is 
seen if the Richardson-Lucy inversion is allow to continue it- 
erating, namely it converges to the small scale noise features 
in the data. Therefore, the smoothest answer is produced 
by the Richardson-Lucy inversion after ~10 iterations. Sec- 
ondly, both of these methods do not do particularly well at 
the end points of the integration range. On the large scales 
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this was to be expected as it is here that the signal is al- 
most zero and the noise dominates. On small scales this 
is due to a combination of the cutoff at 0.1/i~ Mpc and 
the approximation of the integrals in equations [13 and [14 
by finite summations with quite large bins. Thirdly, as a 
consistency check on the Richardson-Lucy method, we have 
examined the Wv(a) produced by the iteration procedure. 
We find that this w„((t) is in excellent agreement with the 
original one obtained directly from the data. Finally, while 
we have only shown the results of testing these inversion 
methods on the CDM mock catalogues, we also tested these 
techniques on other data sets. We constructed simple power 
laws in ^(ct, tt) and also more complicated power laws with 
break features, both with and without noise. The methods 
passed all of these tests and therefore we are confident of ap- 
plying them to the Durham/UKST survey to give believable 
results for £,{r). 

We now discuss the i,{r) results obtained from the 
Durham/UKST survey. Figs, kl and o show that the results 
from these two different inversion techniques agree well. The 
Abel method is generally noisier then the Richardson-Lucy 
method but has ^(r) estimates at twice as many points. 
Comparing the power law fits gives very good agreement in 
both correlation length and slope: ro = 4.8±0.5/i~^Mpc and 
7 = 1.6±0.3 for the Abel method, and ro = 4.6±0.6/i"^Mpc 
and 7 = 1.6 ± 0.1 for the Richardson-Lucy method. These 
values are also in good agreement with those estimated by 
modelling f (r) as a power law and fitting this model to 
U)„((j), namely ro = 5.1 ± 0.3 and slope 7 = 1.6 ± 0.1 
(Paper III). We will use these estimates of ^(r) from the 
Durham/UKST survey in Section ra when comparing the 
real and redshift space correlation functions to estimate a 
value of OP-'^/h. 

To produce a final estimate of the real and redshift 
space 2-point correlation functions we have combined the 
data from the Durham/UKST survey with that of the APM- 
Stromlo survey and the Las Campanas survey (Loveday et 
al. 1992; Loveday et al. 1995; Lin et al. 1996; Tucker et 
al. 1996). These results are shown in Fig. ^. For the red- 
shift space 2-point correlation function (solid circles) we 
have combined the different data sets directly using an error 
weighted mean of the 3 surveys. For the real space 2-point 
correlation function (open circles) we use the Abel inver- 
sion method (for ease) on the 3 different iu„(cr)'s and then 
combine in an error weighted manner. All of the combined 
estimates use the weighting of Efstathiou (1988) and the 
estimator of Hamilton (1993). The error bars assume that 
each survey provides a statistically independent estimate of 
^. For comparison, we plot the real space 2-point correla- 
tion function (solid line) estimated from the inversion of the 
APM w{e) by Baugh (1996). 

We see that the real space ^(r) appears to quite 
nicely approximate a featureless single power law out to 
~20ft~^Mpc. The high point at ~30/i"^Mpc is at the end 
of the inversion region, where this technique becomes unre- 
liable. This estimate agrees well with that of Baugh (1996), 
who also sees a slight 'shoulder' feature on 5 — 25/i~^Mpc 
scales. The redshift space ^(s) appears to have a change of 
shape in both amplitude and slope near ~5/i~^Mpc. The 
smaller scale effect is probably due to the velocity disper- 
sion, which smooths the signal out, while the larger scale 
effect is probably due to the linear infall of galaxies, which 
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Figure 7. Comparison of the real and redshift space 2-point 
correlation functions obtained by combining the data from the 
Durham/UKST survey with that of the APM-Stromlo survey and 
the Las Campanas survey (Loveday et al. 1992; Loveday et al. 
1995; Lin et al. 1996; Tucker et al. 1996). The solid circles show 
the redshift space ^(s), which is an error weighted mean of the 3 
surveys. The open circles show the real space {(r) obtained from 
an error weighted combination of the 3 different Abel inverted 
uii,((t)'s. All of the error bars assume that each survey provides 
a statistically independent estimate of ^. The solid line shows 
Baugh's (1996) real space ^(r) estimate which comes from the 
inversion of the APM w{9). 



amplifies the signal. These effects are discussed in more de- 
tail in Sections n and o. Finally, we make a comment about 
the 'shoulder' feature seen in Baugh's (1996) ^{r). First ev- 
idence for such a feature was observed in ^(s) by Shanks et 
al. (1983) and Shanks et al. (1989) who noted the system- 
atic increase between the amplitude of ^(s) at r ~ 7h~^Mpc 
and the lower amplitudes obtained from Wv{o') at smaller 
scales, as continue to be found here (see Fig. 7). Some of 
this excess in ^(s) is now thought to be due to redshift-space 
distortion from infall (see Section 6) with some contribution 
from a real-space 'shoulder' in ^(r), although the evidence is 
stronger in the 2-D data of Baugh (1996) than in the redshift 
survey data presented here. 



5 NON-LINEAR SCALES: THE 

ONE-DIMENSIONAL PAIRWISE VELOCITY 
DISPERSION 

We now use standard modelling techniques to estimate the 
one-dimensional rms galaxy pairwise velocity dispersion, 
< w^ >^ J, from the Durham/UKST survey. As we saw 
in Section y this velocity dispersion visibly elongates the 
^(ct, tt) contours in the line of sight direction and this is now 
quantitatively measured. 
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5.1 Modelling ^(a,7r) 

We follow the modelling of Peebles (1980) and define v to 
be the peculiar velocity of a galaxy above the Hubble fiow, 
therefore w = Vi — Vj is the peculiar velocity difference 
of two galaxies separated by a vector r. Let g{r, w) be the 
distribution function of w. Therefore 



l + 5(a,^) 



where 



[l + ^{r)]g{r,w)d^w, 



Ho' 



(15) 



(16) 



and Wz is the component of w parallel to the line of sight, 
which for simplicity is called the z direction. It is common 
to assume that g is a slowly varing function of r, such that 
g{r, w) = (?(w). Therefore, one can make the approximation 



(17) 



dwx / dwy g(w) = f{wz)- 
Equation |l5| then becomes 

which further reduces to 
^(cr,7r)= / i(r)f(wz)dwz, 



(18) 



(19) 



when the unit normalisation of f{wz) is considered. A 
streaming model which describes the relative bulk motion 
of galaxies towards (or away from) each other can be incor- 
porated as follows 

g(r,w) =3(w-TO(r)), (20) 

where v(r) is the streaming model in question. Equation M 
then becomes 

1 + C(<7,^)= f [l + i{r)]f{wz-v{rz))dwz, (21) 



f[wz-v{rz)]dwz. (22) 



1 + £. (y^^TTf) 



Obviously, we require models for the real space 2-point cor- 
relation function, £,{r), the distribution function, f{wz), and 
the streaming motion, v{rz). The real space 2-point corre- 
lation function is simply modelled by a power law, ^(r) = 
(ro/r)'', which should be accurate out to ~20/i~^Mpc. For 
the distribution function we tried both exponential and 
Gaussian functions and found that the exponential provided 
a significantly better fit to the shape of the A'^-body simula- 
tion results and so is used exclusively here 



/('".) = 



^/2 < w2 >i/2 



exp 



-^/2. 



< W^ >l/2 



(23) 



where < w'^ >i/2 jg ^j^g j-jjjg pairwise velocity dispersion, 
namely the second moment of the distribution function 



< w > = 



f{wz)wzdwz 



(24) 



A realistic streaming model might be expected to depend 
on the clustering, biasing and mean mass density of the 



universe. The infall model of Bean et al. (1983) takes the 
maximal approach by assuming fl — 1, b = 1 and uses the 
second BBGKY equation (e.g. Peebles 1980) to give 



vivz) 



-Hnr-z 



ar. 



i+ar. 



(25) 



We favour this streaming model in the analysis presented 
here. 



5.2 Testing the Modelling 

We have tested this modelling with both the CDM simu- 
lations and mock catalogues. Before describing these tests 
we first present the answers we are trying to reproduce. We 
have calculated the values of < w^ >i/2 directly from the 



simulations and find that < w > 



1/2. 



950 and 750 kms" 



on ~l/i~^Mpc scales for the SCDM and LCDM models, re- 
spectively. In Paper II we estimated the real space ^(r) cor- 
relation lengths to be ~5.0fe~'^Mpc and ~6.0/i~^Mpc for the 
SCDM and LCDM models, respectively. For both of these 
models the slope of the real space ^(r) was ~2.2. 

In the fitting procedure there are three parameters 
which can be estimated, < w'^ >^", ro and 7. Given our 
available computing constraints we only choose to fit for 
two of these and fix 7 to be a constant value. The results 
obtained from this procedure are insensitive to the value of 
7 chosen, provided a realistic value is used. For the CDM 
simulations we use 7 = 2.2. Also, when the streaming model 
is used we assume that ro = 5.0-6.0/i~^Mpc in equation Ua 
only, again the fits are relatively insensitive to the value 
chosen. We fit < w^ >i/2 ^nd ro using an approximate x^ 
statistic in the range 0-20/i~^Mpc for four different values 
of a, with the standard deviations from the simulations or 
mock catalogues being used accordingly in the x^ statistic. 
The results of these fits using the ^{a, vr) estimated directly 
from the A'^-body simulations are shown in Figs. H and H for 
the SCDM and LCDM models, respectively, with the best 
fit values given in Table uj. The solid histogram denotes the 
averaged ^(cr, tt) in the quoted a range, while the solid and 
dotted lines show the fits with and without the streaming 
model, respectively. Similarly, Table O shows the results of 
the fits to the optimally estimated f (a, tt) from the SCDM 
and LCDM mock catalogues. The quoted error bars on the 
mock catalogue results come from the la standard deviation 
between the mock catalogues and therefore refiects the error 
on an individual mock catalogue. 

The results from the CDM simulations show that 
the streaming model only becomes important (in terms 
of producing consistent results for < w^ >^'^) when 
a > l-2/i~^Mpc. Of course, this assumes that < w^ >i/2 
does not vary with a. The best fit values to the CDM sim- 
ulation results (with streaming), had < w^ >^"= 980 ± 22 
kms'\ ro = 5.00 ± 0.24fe-^Mpc for the SCDM model and 
< w^ >i/2^ g35 _^ 50 kms-\ ro = 5.12 ± 0.69/i"^Mpc for 
the LCDM model. These quoted errors simply come from 
the scatter in the best fit values of Table nl The values of 
can be compared with those estimated directly 



< w > 



1/2 



from the A^-body simulations on ~l/i ^Mpc scales, namely 
950 and 750 kms"^ for the SCDM and LCDM models, re- 
spectively. This agreement is adequate given the slightly 
ad hoc assumption of an exponential distribution function. 
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Figure 8. Histograms of i;{a, it) estimated from the SCDM A'^-body sim ulati ons as a function of vr for different (constant) values of a. 
The solid curve shows the minimum x^ fit using the modelling of Section 5.1 with a streaming model, while the dotted curve shows the 
fit without the streaming model. 



The values of ro can be compared with the approximate real 
space values estimated in Paper III, namely 5.0/i~^Mpc and 
6.0/i"^Mpc for the SCDM and LCDM models, respectively. 
Again, the agreement is adequate in both cases although 
slightly small for the LCDM model. However, closer inspec- 
tion of the LCDM ^(r) plot on r < l/i~^Mpc scales shows 
that ro is slightly lower in this region, ~5.0/i^^Mpc, which 
probably explains our measurement. 



mock catalogue). A similar statement can be made for the 
values of ro estimated from these mock catalogues. 

Overall, our tests confirm that the correct < w^ >i/2 
and ro can be reproduced by this method which models 
the elongation in ^(cr, tt). This holds for both the ^(cr, tt) 
estimated directly from the A'^-body simulations and that 
estimated using the optimal techniques for ^ from the mock 
catalogues, albeit with larger scatter in this case. 



The results from the CDM mock catalogues confirm the 
conclusions drawn from the CDM simulation results. All of 
the < w'^ >i/2 values in Table EI agree well with their coun- 
terparts in Table 111 given the quoted errors (on an individual 



5.3 Results from the Durham/UKST Survey 



We now apply the modelling of Section 5.1 to the 
Durham/UKST survey. We use the optimally estimated 
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Figure 9. Histograms of ^(c, tt) estimated from the LCDM A^-body sim ulati ons as a function of tt for different (eonstant) values of a. 
The solid curve shows the minimum x^ fit using the modelling of Section 5.1 with a streaming model, while the dotted curve shows the 
fit without the streaming model. 



^(a, tt) dat a which was shown in Fig. H and described in 
Section 3.3. Fig. hol shows the histograms of ^(cr, tt) as a func- 
tion of TT for different (constant) values of a. The so lid c urve 

with 



5.1 



shows the minimum x fit of the model from Section 
the streaming model. The dotted curve shows the same fit 
without the streaming model. The results of these fits are 
given in Table pi where we have assumed 7 = 1.7. Small 
changes in this value of 7 do not significantly affect our 
results. We use the standard deviations on an individual 
LCDM mock catalogue in the x^ fits. These fits are more 
than likely biased low by the non-independent nature of the 
points. The errors on the parameters quoted in Table H come 
from the appropriate Ax^ contour about the minimum. 



We can make a few comments about these results. 
Firstly, the streaming model is again required to produce 
consistent fits for < w^ >i/2 fQj- ^ > l-2/i~^Mpc (assum- 
ing that it does not change with a). Secondly, although not 
shown we also fitted for a Gaussian distribution function 
(rather than the exponential one favoured by the simula- 
tions) and found that the noise in the data implied a sim- 
ilar quality of fit for both forms. Therefore, this data did 
not allow one to distinguish between the two models. Fi- 
nally, we combine the results for the exponential velocity 
distribution function with the streaming model at each a 
value to produce an overall estimate of < ui^ >i/2 ^nd ro- 
Assuming that each quoted estimate is independent gives 
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Table 1. Minimum x^ fits for < w^ >l/2 and ro from the CDM 
Af-body simulations using the exponential distribution function 
with and without the streaming model. 



a 
(/i-iMpc) 



< W" > 



1/2 



(kms ^ ) 



(/i^iMpc) 



SCDM 
Streaming Model 



[0,0.5] 


990 




5.1 


[0.5,1] 


1000 




5.3 


[1,2] 


980 




4.8 


[2,4] 


950 




4.8 




No St] 


reaming 


IVIodel 


[0,0.5] 


970 




5.1 


[0.5,1] 


960 




5.4 


[1,2] 


800 




4.9 


[2,4] 


550 


LCDM 


5.1 




Streaming JVIodel 


[0,0.5] 


770 




4.1 


[0.5,1] 


810 




5.3 


[1,2] 


850 




5.4 


[2,4] 


910 




5.6 




No St] 


reaming 


Model 


[0,0.5] 


780 




4.2 


[0.5,1] 


760 




5.4 


[1,2] 


720 




5.6 


[2,4] 


570 




5.9 



X 



(Nfc,;„ = 40) 



40.0 
24.7 
21.2 
34.4 

46.3 
32.3 

28.7 
48.2 



13.8 
7.8 
16.5 
11.1 

14.2 
8.1 
14.0 
9.0 



Table 3. Minimum x^ fits for < w^ >i/2 and ro from the 
Durham/UKST survey using the exponential distribution func- 
tion with and without the streaming model. 



a 


< w'^ >^/^ ro 


x' 


(/i-iMpc) 


(kms-i) (h-^Mpc) 

Durham/UKST 
Streaming Model 


{^Mn = 


[0,0.5] 


510 ±120 5.1 ±0.6 


19.0 


[0.5,1] 


300 ±50 4.7 ±0.3 


23.5 


[1,2] 


500 ±65 4.5 ±0.2 


29.6 


[2,4] 


500 ±65 4.7 ±0.4 
No Streaming Model 


48.0 


[0,0.5] 


470 ±130 5.2 ±0.6 


19.1 


[0.5,1] 


180 ±70 4.8 ±0.4 


24.0 


[1,2] 


270 ±90 4.8 ±0.3 


29.5 


[2,41 


180 ±80 5.5 ±0.2 


56.1 



40) 



< w^ >i/2 



416 ± 36 kms-\ ro = 4.6 ± 0.2/i-^Mpc. We 
make the comment that this value of ro agrees well with 
those estimated from the Durham/UKST survey using other 
methods, see Section nA. These results are discussed in more 
detail in Section M where comparisons with the results from 
structure formation models and other redshift surveys are 
made. 



Table 2. Minimum x^ fits for < w^ >l/2 and ro from the CDM 
mock catalogues using the exponential distribution function with 
and without the streaming model. 






ro 



(h^^Mpc) (kms 



(h-^Mpc) 



SCDM 

Streaming Model 

[0,0.5] 930 ±190 5.0 ±0.4 

[0.5, 1] 970 ± 180 5.3 ± 0.4 

[1,2] 920 ±260 4.6 ±0.5 

[2,4] 810 ±180 4.5 ±0.6 

No Streaming Model 

[0,0.5] 920 ±190 5.0 ±0.4 

[0.5,1] 910 ±170 5.3 ±0.4 

[1,2] 750 ±230 4.8 ±0.4 

[2,4] 450 ±130 5.0 ±0.4 

LCDM 

Streaming Model 

[0,0.5] 610 ±180 3.8 ±0.3 

[0.5,1] 790 ±230 5.2 ±0.7 

[1,2] 840 ±130 5.2 ±0.7 

[2,4] 780 ±150 5.2 ±0.7 

No Streaming Model 

[0,0.5] 590 ±170 3.9 ±0.3 

[0.5, 1] 750 ± 220 5.2 ± 0.6 

[1,2] 690 ±130 5.3 ±0.7 

[2,4] 440 ±140 5.6 ±0.6 



6 LINEAR SCALES: INFALL AND Q.°-^/B 

We now use modelling techniques developed from the lin- 
ear theory result of Kaiser (1987) to estimate the quantity 
Q'^'^ /b from the Durham/UKST survey. We saw in Section H 
that on non-linear scales the velocity dispersion was mainly 
responsible for the anisotropies in ^(cr, tt). However, in Sec- 
tion H we saw that to produce consistent results it was also 
necessary to incorporate a model which imitated the stream- 
ing motions of galaxies. It is this infall of galaxies into over- 
dense regions that visibly distorts the ^{a, n) contours on 
larger scales and this compression in the line of sight direc- 
tion is now quantitatively measured. 



6.1 Modelling ^(a,7r) 

Kaiser (1987) used the distant observer approximation in 
the linear regime of gravitational instability to show that 
the strength of an individual plane wave as measured in 
redshift space is amplified over that measured in real space 
by a factor 

K = 5l{i + ^^lJ{n)/b) , (26) 

where 51- and 5^ are the Fourier amplitudes in real (r) 
and redshift (s) space, respectively, /Xki is the cosine of 
the angle between the wavevector (k) and the line of sight 
(1), /(n) — O,^'^ is the logarithmic derivative of the fiuc- 
tuation growth rate (e.g. Peebles 1980) and b is the lin- 
ear bias factor relating the galaxy and mass distributions, 
(Ap/p)g = b{Ap/p)m- The distant observer approximation 
restricts use of equation M to opening angles in a redshift 
survey of less than ~50° which can cause a systematic effect 
at the ~5% level in f{n)/b (Cole, Fisher & Weinberg 1994). 



© 0000 RAS, MNRAS 000, 000-000 



The Durham/ UKST Galaxy Redshift Survey 13 




10 

TT (h^'Mpc) 

Figure 10. Histograms of ^(cr, n) estimated from the Durham /UKST sur vey as a funetion of tt for different (constant) values of a. The 
solid curve shows the minimum x^ fit using the modelling of Section 5.1 with a streaming model, while the dotted curve shows the fit 
without the streaming model. 



This (1 + /i^i/(r2)/6) factor propagates directly through to 
the power spectrum, P{k,fj,k\) = (SkS^.) 



P'(fc,Mki) = P'Xfc)(l + Mki/(«)/6)' 



(27) 



where the real space P^ik) is assumed to be an isotropic 
function of k only. Thus the anisotropy is a strong function 
of the angle between k and 1. Integrating over all /iki gives 
the angle- averaged redshift space P^ik) 



P'ik) 



/_^jd/iki-P°(fe,Mki) 
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(28) 



= P'-(fc) 



3 6 5 \ b 



(29) 



Fourier transforming equation ES (which has no explicit /iki 
dependence) gives the corresponding relation between the 
angle- averaged ^(s) and ^(r) 



C(s) = CM 



2/(») 1 //(») 
3 b 5\ b 



(30) 



assuming that ^(r) is an isotropic function of r only. If the 
volume integral of ^ is defined as 
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J'i{x) = / £.(y)y dy, 

Jo 
then 



(31) 



J3(s) = Jair) 



1 + 



2 /(n) 

3 b 



l(^) 



(32) 
as two of our methods for 



We will use equations pfl and 
estimating r2°'^/b. 

Hamilton (1992) has extended Kaiser's (1987) analysis 
by Fourier transforming equation M and noting that the 
cosine factor in Fourier space, /i^i = kf/k'^, becomes a dif- 
ferential operator in real space, (9/9jlj)^(V^)^^. Therefore 
equation ET^ becomes 

e{r,^.,) = {l + {f{n)/b){d/d\l\)\\7Y'ye{r), (33) 

where /iri is the cosine of the angle between the pair separa- 
tion, r, and the line of sight, 1. Hamilton (1992) then shows 
that the solution of this equation can be written in terms 
of the first three even spherical harmonics of ^°(r, /iri), with 
the higher moments zero (and all odd moments zero by def- 
inition) 

r(nMrl) = ?o(r)Po(Mrl) +6(r)P2(Mrl) +C4(r)P4(/irl), (34) 

where Pi{jj,ri) are the usual Legendre polynomials and 

2; + i 



^lir) = 



(.''ir,flrl)Pl{fJ'rl)dflr 



(35) 



Arguably the most useful form of Hamilton's (1992) solution 
is that expressed in terms of ^o and ^2 only 



1 + 



2 /(n) 



H^Y 



or by defining 

io{r) = -Co(0 + ^ 

|2(r) = -C2(r), 
it can be written as 

Co 



3 

73- 



i;; ^o{s)s^ds) 



^o(s)s ds, 



[1 + Il3+If3^)' 



(36) 

(37) 

(38) 



We will use the ratio in equation pq as our third method for 
estimating OP'^' /b. 



6.2 Testing the Modelling 

We have tested this modelling with the LCDM simulations 
and mock catalogues. The SCDM model was not analysed 
because we felt that the one-dimensional rms galaxy pair- 
wise velocity dispersion too strongly dominated the £,{a, tv) 
plots for a sensible answer to be obtained. As will be shown 
below this is also the c ase for some aspects of the LCDM 
model. Also, in Section 6.1 we noted that equation Eq and 



the subsequent analysis is only strictly correct in the distant 
observer approximation and that data from redshift survey 
opening angles of < 50° should only really be used. For red- 
shift surveys geometrically similar to the Durham/UKST 
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Figure 11. Estimates of U'-^-^ /b obtained from the LCDM model 
using the method involving the ratio of real to redshift space 2- 
point correlation functions (equation pQ). The solid line denotes 
the results obtained by averaging the ^'s from each A'^-body sim- 
ulation and then manipulating to get an estimate of Q^'^ /b. The 
shaded area denotes the 68% confidence region in Sl^'^/fe obtained 
from the scatter seen between the TV-body simulation ^ results. 
The solid points denote the mean of the results obtained from the 
mock catalogues. The error bars on these points are the observed 
standard deviation on an individual mock catalogue. 



survey this restriction has a negligible impact on the analy- 
sis techniques. 

In Figs, hll, h2 and hs the following conventions are 
adopted: the dotted line denotes the expected value of 
fin)/b ~ n°-^/b = (0.2)°Vl - 0.38; the solid line denotes 
the result for fi^'^/b obtained fronr the average of the LCDM 
simulations (i.e. estimate an averaged ^(cr, tt) from the 5 sim- 
ulations and then manipulate to get one value of OP'^ /b); the 
shaded area denotes the 68% confidence region in ff'^/b ob- 
tained from the scatter seen between the LCDM simulations 
(i.e. take the ^{a, n) from each of the 5 simulations, manip- 
ulate to get 5 values of Q}''^ /b and then average at the end); 
and the points with error bars are the mean and la scat- 
ter seen in the LCDM mock catalogues (i.e. take the ^(cr, tt) 
from each of the 15 mock catalogues, manipulate to get 15 
values of Q''^'^' /b and then average at the end). These errors 
are the observed standard deviation between the mock cata- 
logues and therefore reflect the error on an individual mock 
catalogue. 

The results from equation BG are shown in Fig. [01. This 
method uses the ratio of the real to redshift space ^'s to 
estimate ff'^/b. We estimate ^(s) using the methods de- 
scribed in Paper HI. We estimate ^(r) using the Abel in- 
version of the projected correlation function, iu„ (a), with 
a TTcut = 30/i~^Mpc, as described in Section ^ Examining 
the results obtained from the averaged ^ of the simulations 
(solid line), it appears that the method itself is not domi- 
nated by non-linear effects above ~6/i~^Mpc, although they 
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Figure 12. The same as Fig. |ll| but the Q'^-^/6 results use the 
method involving the ratio of real to redshift space volume inte- 
grated 2-point correlation functions (equation B2|) . 



sensitive to non-linear effects but the large scatter (for the 
mock catalogues) renders the method almost useless in this 
case. 

The results from equation p3 are shown in Fig. [13. This 
method uses the ratio of the real to redshift space Ja's to 
estimate n^'^/fc. We estimate Jsis) by evaluating the inte- 
gral in equation pl| using the above ^{s). We estimate Jair) 
from the same integral but using the above f (r). Examining 
the results obtained from the averaged ^ of the simulations, 
it appears that the method itself is not dominated by non- 
linear effects above ~15/i~^Mpc. The results obtained by 
averaging the 5 estimates of ff'^/b from each simulation 
are similarly consistent with the expected value. The results 
obtained from the mock catalogues also reproduce the cor- 
rect answer but with a larger scatter. Overall, this method is 
considerably more successful (in terms of more accurate re- 
sults) than the (,{s)/£,{r) one. However, one does pay a price 
by having to go to larger scales before the non-linearities 
become negligible and the integration procedure also makes 
the points highly correlated. 

The results from equation pq are shown in Fig. uA This 
method uses the ratio involving the second and zeroth spher- 
ical harmonics off to estimate OP'^ /b. The ^;'s are estimated 
using equation Ba, which in practice becomes 



^iir) = {2l + l)Afiri 



{Mi,>0) 

E 



^%r,^li)Pi{fj.l 



(39) 
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Figure 13. The same as Fig. nil but the fl'^-^/b results use the 
method involving a ratio incorporating the second and zeroth 
spherical harmonics of the redshift space 2-point correlation func- 
tion (equation 38h. 



could cause the ~0.1 systematic offset in Q,^'^ /b that is seen 
out to > 30/i~^Mpc. Unfortunately, the results obtained by 
averaging the 5 estimates of Q°'^ /b from each (full) simula- 
tion show that noise begins to dominate the inversion pro- 
cess between 15-20/i~^Mpc. For the mock catalogues noise 
dominates at all scales and the results resemble a scatter 
plot in places. Overall, this method appears relatively in- 



and we use a binning of A^ri = 0.2. Examining the re- 
sults obtained from the averaged ^ of the simulations, it 
appears that the method itself is severely affected by non- 
linearities which cause a negative value of OP'^ /b to be mea- 
sured until ~13/i~^Mpc. (Such a negative value is obviously 
unphysical and simply due to the shape of the ^(cr, tt) con- 
tours.) Similar results are found by averaging the 5 esti- 
mates of Q,"'^ /b from each simulation. The results obtained 
from the mock catalogues trace the simulation results ad- 
equately, although not on scales < 10/i~^Mpc where they 
appear biased high. Overall, the impression is that the non- 
linearities make this method redundant. However, when con- 
sidering the Durham/UKST survey one must remember that 
its value of the velocity dispersion is almost half that of the 
LCDM model. Therefore, the elongation in ^(cr, tt) which 
causes the above problems will be significantly lower. 



6.3 Results from the Durham/UKST Survey 



We now apply the modelling of Section 3.1 to the 



Durham/UKST survey. We calculate our results using the 
optimally estimated 2-point correlation function described 
in Section 3.3 and Paper HI. Before presenting our results 
of fi^'^/fo we show the data which produced them. Fig. u4 
shows the real and redshift space ^'s, where £,{s) is estimated 
directly from the redshift survey and ^(r) comes from the 
Abel inversion of Section [4.2| shown in Fig. 0. Fig. uM shows 
the real and redshift space Js's, where the integral in equa- 
tion Kll is applied to the above £,'s. Fig. uM shows the zeroth 
and second spherical harmonics estimated using equation ^ 
on the optimally estimated ^''(r, pri)- The error bars shown 
are those estimated from the scatter between the 4 quad- 
rants of the Durham/UKST survey. These errors were of a 
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Figure 14. Comparison of the real and redshift space 2-point cor- 
relation functions optimally estimated from the Durham/UKST 
survey, ^(s) was estimated directly from the redshift survey (solid 
triangles connected by th e do tted line) and ffr) comes from the 
Abel inversion of Section [4. 2] shown in Fig. W (solid circles con- 
nected by the solid line). The error bars denote the scatter seen 
between in the 4 quadrants of the Durham/UKST survey. 



very similar size to those estimated on an individual LCDM 
mock catalogue. 

We comment briefly on the results in each of these 
figures. Firstly, Fig. h4^ shows that ^(r) > ^(s) below 
~l/i~^Mpc, while g( s) > ^(r) on larger scales. This was 
Unfortunately, the noise in these 



described in Section 4.4 



estimates is very high and it will more than likely make any 
measurement of OP'^ /b redundant. Secondly, Fig. ha shows 
that, once again, at small separations Jd.ir) > J3(s), while 
at larger separations JsCs) > Js(r-). There appears to be 
a near constant offset in the real/redshift Ig Js's on scales 
10-20/i~^Mpc (i.e. a constant multiplicative factor in lin- 
ear J3). This should give a consistent estimate of OP'^ /b on 
these scales. Finally, Fig. Ilq shows that ^2 is positive until 
~8/i~'^Mpc, which is caused by the elongation of the ^(cr, tt) 
contours parallel to the line of sight. On larger separations 
^2 is negative due to the compression of the (^{a, tt) contours 
parallel to the line of sight. Therefore, there appears to be 
a signiflcant signal to measure in this method. 

We now present the results for OP'^ /b estimated from 
the Durham/UKST survey. Fig. hj shows the results of 
applying equation B(j (solid squares), equation B2^ (solid 
triangles) and equation BS (solid circles) to the data in 
Figs. [I4I, [la and |16l respectively. The plotted errors de- 
note the standard deviation seen between the 4 quadrants 
of the Durham/UKST survey. Similarly sized errors were 
seen in the scatter between the individual LCDM mock cata- 
logues. For clarity, error bars are not shown for the (,{s)/^{t) 
method because they are very large and only cause confu- 
sion. 

We can make a few comments about the results from 
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Figure 15. Comparison of the real and redshift space Js's es- 
timated from the Durham/UKST survey. These were calculated 
using the integral in equation Kll with the ^'s of Fig. [14. Again, 
the error bars denote the scatter seen between in the 4 quadrants 
of the Durham/UKST survey. 
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Figure 16. The zeroth and second spherical harmonics from t he 
Durham/UKST survey. These were calculated using equation p9 
on the optimally estimated ^''(r, /^ri)- The solid circles connected 
by the solid line shows ^Oi while the solid triangles connected by 
the dotted line shows —li,2- Note that — §2 is plotted so that it 
is positive in the large scale region of interest. Again, the error 
bars denote the scatter seen between in the 4 quadrants of the 
Durham/UKST survey. 
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Figure 17. Estimates of Q^-^/b obtained from the Durham/ 
UKST survey. We use the three methods from equation BCJ (soUd 
squares), equation p3 (solid triangles) and equation pq (s olid cir- 
cles) and apply them to the data in Figs. h4|, lla andTa, respec- 
tively. The plotted errors denote the standard deviation seen be- 
tween the 4 quadrants of the Durham/UKST survey. For clarity, 
error bars are not shown for the ^(s)/^(r) method because they 
are very large and only cause confusion. 



these three methods. Section |6.2 
interest is 



showed that our region of 
10-30/i~^Mpc because of non-linear effects on 
smaller scales and noise on larger scales. For the (,{s)/^{r) 
method, the points have no systematic trend (other than a 
large random scatter) and we discount them from the fur- 
ther analysis. For the J3(s)/J3(r) method, we obtain quite 
consistent results although one must remember that these 
points are not independent because of the integration prod- 
edure. We choose to quote the value at ~20/i~^Mpc as being 
representative of our results, Q*^'^ /b — 0.52 ± 0.39. We be- 
lieve that we are being very conservative with this quoted 
error and, for example, the point at ~12/i~^Mpc only has an 
error of ±0.19. For the C2/C0 method, we again obtain very 
consistent results and quote the value at ~18/i^^Mpc as be- 
ing representative, OP'^ /b = 0.45 ± 0.38, where the quoted 
error is an average of the error bars in the ~10-30ft~'^Mpc 
region and is therefore the typical error on any individual 
point in this region. Naively combining the results in a max- 
imum likelihood manner in the ~10-20/i^^Mpc region gives 
an estimate of Q''^'^' /b = 0.48 ± 0.11. In this case the error 
estimate is likely to be slightly underestimated given the 
generally non-independent nature of ^ points. These results 
are discussed in more detail in Section M where comparisons 
with the results from structure formation models and other 
redshift surveys are made. 



7 DISCUSSION 

We now discuss the results obtained from the analysis pre- 
sented in Sections H and M. Given the observed problems with 
the unweighted estimate of ^ (Paper III), we now favour the 
weighted estimate of ^ in all of our analysis. This explains 
any slight numerical differences with respect to the quoted 
results of Ratcliffe et al. (1996a). However, we state clearly 
that our overall conclusions remain completely unchanged. 



7.1 The One-Dimensional RMS Pairwise Velocity 
Dispersion 

The minimum x^ fit of an exponential distribution func- 
tion (with a streaming model) to the optimally estimated 
(,{a, n) from the Durham/UKST survey gave a value for the 
one-dimensional rms pairwise velocity dispersion of 416 ±36 
kms~^. We note that this quoted error is the formal er- 
ror obtained from the x^ statistic. Therefore, is it likely 
to be slightly underestimated given the non-independent 
nature of ^. Our observed value is of particular interest 
as recent estimates from new redshift surveys and the re- 
analysis of old redshift surveys have been measuring larger 
velocity dispersions than the canonical value of 340 ± 40 
kms"'^ found by Davis & Peebles (1983) from the CfAl sur- 
vey. For example, using the CfA2/SSRS2 survey, Marzke et 
al. (1995) find 540 ± 180 kms~^ and, using the Las Cam- 
panas survey, Lin et al. (1996) find 452 ± 60 kms^^. Also, 
Mo, Jing & Borner (1993) measured large variations (200- 
1000 kms"'^) in the velocity dispersion for a number of sam- 
ples of similar size to CfAl and they show that the esti- 
mated velocity dispersion is sensitive to galaxy sampling, 
especially dominant clusters the size of Coma. This new 
estimate from the Durham/UKST survey is still on the 
slightly low side, supporting the old Davis & Peebles (1983) 
value, but is not inconsistent (>3(t) with any of these other 
measured values. When considering these values it is im- 
portant to note that the Durham/UKST survey covers a 
volume ~4 x 10^h~'^Mpc'^, approximately twice that of the 
CfA2/SSRS2 survey and half that of the Las Campanas sur- 
vey. Also, in an unbiased (COBE-normalised) CDM model, 
Marzke et al. (1995) estimated that the velocity dispersion 
would converge to 10% within a volume ~5 x 10'°h~^Mpc'^. 
Therefore, the measurement from the Durham/UKST sur- 
vey is hopefully both believable and representative of the 
actual value in the Universe. Finally, we note that the 
Durham/UKST survey does not contain any extremely dom- 
inant clusters (of Coma-like size) and therefore will not be 
biased high by this. 

The best estimate of the one-dimensional rms pairwise 
velocity dispersion from the SCDM and LCDM simulations 
was 980 and 835 kms~^, respectively, obtained using an ex- 
ponential distribution function. These estimates agree well 
with the actual value of the velocity dispersion measured di- 
rectly from the A'^-body simulations. However, both of these 
values are inconsistent with the measured value from the 
Durham/UKST survey at high levels of significance. In fact, 
taking the most negative approach possible (i.e. using the in- 
dividual mock catalogue velocity dispersion error bar), one 
still finds a significant rejection (>3a) of both CDM mod- 
els. However, it should be noted that a significant velocity 
bias, bv, between the matter and galaxy velocity distribu- 
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tions (6^~0.4), see Couchman & Carlberg (1992), would 
allow consistent results between the models and the data. 
Also, this rejection of the CDM models assumes that the 
simple models of linear biasing used here (Bardeen et al. 
1986) are an adequate description of the galaxy formation 
process. 

7.2 Infall and n°-^/b 

The best estimates of ri""/b from the Durham/UKST sur- 
vey are 0.52±0.39 for the J3(s)/J3(r) method and 0.45±0.38 
for the ^2 /Co method, where we have (conservatively) quoted 
the estimate at one separation only. Given the integration 
procedure involved in the J3 method we do not attempt to 
combine these points. However, naively combining the C2/C0 
estimates gives OP'^ /b = 0.48 ± 0.11. This error is likely to 
be slightly underestimated given the non-independent na- 
ture of 5. Our estimates can be compared with other op- 
tical values of Vp'^ /h estimated using similar methods in- 
volving redshift space distortions. Peacock & Dodds (1994) 
use the real and redshift space power spectrum estimates 
of various cluster, radio, optical and IRAS samples to mea- 
sure n"*^/?) = 0.77 ± 0.16. Loveday et al. (1996) use the 
Ja method to measure n"'^/fe = 0.48 ± 0.12 for the APM- 
Stromlo survey. Lin et al. (1996) use the ^2/^0 method to 
measure 0°®/& = 0.5 ± 0.25 for the Las Campanas survey. 

As one can see, all these observed values of OP'^ /b are 
consistent with ~0.5. Therefore, using two fiducial values 
for h, b = 1 implies and the universe is open, f2~0.3, and 
6~2 implies the universe has the critical-density, f2 = 1. 

Finally, we state that our OP'^ /b estimates are consis- 
tent with those from the two CDM models of structure for- 
mation considered here (SCDM and LCDM). However, given 
that these models predict (f2'''®/6)~0.4-0.6 we are unable to 
distinguish between them. 



8 CONCLUSIONS 

We have investigated the redshift space distortions in the 
Durham/UKST Galaxy Redshift Survey using the 2-point 
correlation function perpendicular and parallel to the line 
of sight, C(o-, 7r). On small, non-linear scales we observe an 
elongation of the ^ contours in the line of sight direction, 
which is due to the velocity dispersion of galaxies in virial 
regions. This is the common "Finger of God" effect seen 
in redshift surveys. On larger, linear scales we observe a 
compression of the ^ contours in the line of sight direction, 
which is due to the infall of galaxies into overdense regions. 
We attempt to estimate the real space 2-point corre- 
lation function by direct inversion of the projected corre- 
lation function. We use two difi'erent inversion methods, 
Abel inversion and a new application of the Richardson- 
Lucy technique. We have tested these methods on mock 
catalogues drawn from cold dark matter (CDM) A^-body 
simulations and find that they reproduce the correct an- 
swer. We apply the methods to the Durham/UKST sur- 
vey and estimate £,{r). We find that a simple power law 
model gives best fit real space parameters of the correla- 
tion length. To — 4.8 ± 0.5fe~^Mpc for the Abel method and 
ro — 4.6±0.3/i^^Mpc for the Richardson-Lucy method, and 
slope, 7 = 1.6 ± 0.3 for the Abel method and 7 = 1.6 ± 0.1 



for the Richardson-Lucy method. Our estimate is consistent 
with those from other redshift surveys and with the inver- 
sion of the APM w{e) by Baugh (1996). 

We use standard modelling techniques (e.g. Peebles 
1980) to estimate the one- dimensional rms pairwise veloc- 
ity dispersion of galaxies in the Durham/UKST survey. 
These methods (which were again tested on the mock cat- 
alogues) give an estimate of < w'^ >i/2_ 4^5 -|- 35 kms~^ 
on ~l/i^^Mpc scales. This value agrees well with recent es- 
timates from new redshift surveys and re-analysis of old 
redshift surveys (Markze et al. 1995; Lin et al. 1996; Mo 
et al. 1993), although our value is still consistent with the 
canonical value of Davis & Peebles (1983). We compare 
with the predictions of the standard CDM model (SCDM; 
Qh — 0.5 &i b = 1.6) and the low density CDM model with 
a non-zero cosmological constant to ensure spatial flatness 
(LCDM; D,h = 0.2, A = 0.8 & 6 = 1.0). We find that our re- 
sults are significantly (>3cr) below the estimates from these 
models (assuming that linear biasing applies) and therefore 
these models are inconsistent with the data in this con- 
text. This modelling also produces an estimate of the real 
space correlation length from the Durham/UKST survey of 
ro = 4.6 ± 0.2/i~^Mpc, which is very consistent with the 
previously quoted values. 

We use the modelling techniques of Kaiser (1987) 
and Hamilton (1992) to estimate OP'^/b from the 
Durham/UKST survey. We test the methods with the 
mock catalogues and find that consistent results can be 
obtained. The Durham/UKST survey results from the 
(,{s)/(,{r) method are too noisy to obtain a useful answer. 
The Js{s)/J;i{r) method produces an estimate of OP'^/b — 
0.52 ± 0.39, where we quote the estimate (and error) at one 
point only (~20ft^^Mpc) due to the non-independent na- 
ture of the integration procedure. The C2/C0 method gives 
n""/?) = 0.45±0.38 at ~18/i"^Mpc, which is representative 
of the results on these scales. Naively combining the points 
in the range ~10-20/i~^Mpc gives a maximum likelihood fit 
of 0.48 ± 0.11, where the error is likely to be slightly under- 
estimated due to the general correlated nature of ^ points. 
A comparison with other optical estimates of OP'^ /b from 
redshift space distortion methods gives very consistent re- 
sults with a value of (r2'''®/6)~0.5. This argues against an 
unbiased critical-density universe {b = 1 &i fl — 1), instead 
favouring either an unbiased low density universe {b — 1 &i 
f2~0.3) or a biased critical-density universe (fe~2 & fi = 1). 
Also, given that both CDM models considered here predict 
(n°''/fe)~0.4-0.6 we cannot use our r2°'^/fo resuhs to distin- 
guish between them. 

Overall, combining these results with those presented in 
Paper III, we find that the standard CDM model underpre- 
dicts our 2-point correlation function results at large scales 
and overpredicts the one-dimensional pairwise velocity dis- 
persion at small scales. Therefore, our results argue for a 
model with a density perturbation spectrum more skewed 
towards large scales, such as a low n CDM model with a 
cosmological constant. 
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